Microscopic Superfluidity in Bose Gases: From 3D to ID 
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The superfluid fraction of ideal and interacting inhomogeneous Bose gases with varying asymmetry 
is investigated at finite temperature using well-known properties of the harmonic oscillator as well 
as the essentially exact microscopic path integral Monte Carlo method. We find that the superfluid 
fraction (i) is essentially independent of the interaction strength for all temperatures considered, 
(ii) changes profoundly as the effective dimensionality is varied from three- to one- dimensional, (iii) 
is approximately equal to the condensate fraction No/N for spherical Bose gases, and (iv) deviates 
dramatically from No/N for highly-elongated Bose gases. 



Macroscopic objects such as liquid 4 He show many pe- 
culiar properties that can be attributed to superfluid- 
ity Q. Among these are the absence of viscosity, the 
occurence of persistent currents, the existence of vortices 
and the reduction of the moment of inertia. Connec- 
tions between manifestations of superfluidity and Bose 
Einstein condensation have been studied extensively in 
the context of liquid 4 He since the discovery of its super- 
fluidity in 1938. While much progress has been made in 
our understanding of such strongly interacting systems, 
many questions remain unanswered. 

Thanks to the realization of gaseous Bose Einstein con- 
densates in 1995 0, the study of superfluid effects of 
mesoscopic systems has become possible. Indeed, the 
creation of vortices Q and vortex lattices has been 
demonstrated in inhomogeneous Bose gases, and, most 
recently, also in degenerate Fermi gases in the BEC-BCS 
crossover regime |5j. Following the work on 4 He enclosed 
in a cylinder , superfluidity of inhomogeneous systems 
can be described through their rotational properties. The 
superfluid fraction is defined by the departure of the 
quantum mechanical moment of inertia 0^ with respect 
to h from its classical, or rigid, value 6" ff . Here, the 
moment of inertia O^, @n = (d(L ■ n) u /du)) u= Q, is de- 
fined by the linear response of the system to a rotational 
field H ext — —uj ■ L, where uj = ton; to denotes the angu- 
lar frequency and L the total angular momentum. The 
thermal expectation value (-^ is evaluated for the sys- 
tem perturbed by H ext . The normal fraction is the part 
of the system that responds classically, i.e., Qa/Q^ 9 , and 
the superfluid fraction (n s /n)n is 1 — Q^/Q 1 ? 3 . 

The temperature dependence of Qfi has been evaluated 
for a non-interacting Bose gas under harmonic confine- 
ment in the so-called macroscopic approximation . The 
effects of weak interactions have been estimated within 
the Thomas-Fermi approximation Q ■ This Letter deter- 
mines the temperature-dependence of 0^ and O^ 9 , and 
hence of the superfluid fraction, for small atomic gases 
with N = 27 bosons for varying confinement and inter- 
action strength non-perturbatively using the essentially 
exact microscopic path integral Monte Carlo (PIMC) 
method ||. For the ideal gas, we additionally deter- 
mine thermal expectation values using well-known prop- 



erties of the harmonic oscillator M. In contrast to 4 He 
clusters ^3 or deformed nuclei ^lj, whose interaction 
strength and internal temperature are largely "set by na- 
ture" , atomic gases provide us with unprecedented con- 
trol. The temperature can be controlled by changing the 
cooling scheme 12] , the interaction strength can be tuned 
by applying an external mag netic field in the vicinity of a 
Feshbach resonance Il4j |. and the dimensionality can 
be reduced by varying the external confinement |l5l llfi| . 
Here, we focus on the crossover from three-dimensional 
(3D) to one-dimensional (ID) behavior. We show that 
reduced dimensionality leads to an increase of the su- 
perfluid response. The superfluid fraction for 3D gases 
roughly coincides with the condensate fraction Nq/N. 
In the quasi-lD regime, however, the superfluid fraction 
is much larger than No/N. Our calculations show that 
the superfluid response depends, if at all, weakly on the 
strength of the atom-atom interactions. 

Consider N bosons with mass m under external har- 
monic confinement, 
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Here, pj and Zj denote the transverse and longitudinal 
coordinate of the jth atom, respectively, and oj p and lu z 
the transverse and longitudinal angular frequency of the 
trapping potential, respectively. The atom-atom poten- 
tial V depends on the interparticle distance rjk between 
atom j and atom k. For the non-interacting gas, i.e., 
V(r) — 0, we calculate thermal expectation values in the 
grandcanonical ensemble using well-known properties of 
the harmonic oscillator |9j. To simulate effectively re- 
pulsive Bose gases, we use a hard sphere potential V(r) 
with 3D atom-atom scattering length a; in particular, 
a = 0.00433 and 0.0433a z , where a z — \Jh/ (mu z ). In 
this case, we use the numerically more involved PIMC 
technique , which determines thermal expectation val- 
ues in the canonical ensemble. For the purpose of the 
present study, differences between expectation values cal- 
culated in the grandcanonical and in the canonical en- 
semble (see also Ref. 0) are negligible. 

To investigate the crossover from 3D to ID for N 
bosons, we vary the angular frequency uj p such that 



2 




20 40 60 80 100 
Aspect ratio L 



FIG. 1: Approximate 3D transition temperature fcsT c (see, 
e.g., Eq. (19) of Ref. [T^) in units of hu z (ui z — 2irv z ) as 
a function of L for TV = 27. Vertical arrows indicate the 
interval 0.1 < T/T c < 1.4 for L = 1, 10 and 100. The inset 
shows the specific heat C divided by NUb calculated in the 
grandcanonical ensemble as a function of T/T c for TV = 27 
and L = 1 (solid line), 10 (dotted line) and 100 (dashed line). 



L = 1, 10 and 100, where L = lo p /uj z . The approxi- 
mate 3D transition temperature T c , obtained for vanish- 
ing atom-atom interactions, then depends on u> p , io z and 
N (the T c used throughout this paper includes finite-size 
corrections; see, e.g., Eq. (19) of Ref. |l8|), A dotted line 
in Fig. ^ shows T c as a function of the aspect ratio L for 
TV = 27. Below, we report calculations for L = 1, 10 and 
100 over a wide temperature range, i.e., 0.1 < T/T c < 1.4 
(see vertical arrows in Fig.^). 

Highly-elongated gases at T = can to a very good 
approximation be described by an effective ID Hamil- 
tonian for any 3D scattering length a if N/L <C 1 [l9j . 
For TV = 27 and L = 100, we find N/L = 0.27. At fi- 
nite temperature, the behavior of highly-elongated Bosc 
gases depends on two energy scales, the oscillator en- 
ergies hujp and hu) z of the tight and weak confinement 
direction, respectively. For N = 27 and L — 100, three 
temperature regimes exist [2(| : i) T is larger than the 3D 
transition temperature T c (excited transverse modes are 
occupied); ii) T is lower than T r but larger than the ID 
transition temperature T C 1D 21] (transverse excitations 
are largely frozen out); and iii) T is smaller than T c 1£) 
(excited longitudinal modes are largely frozen out). For 
N = 27 and L = 100, the approximate 3D transition 
temperature is fe^Te = 36.0?kt; z , while the approximate 
ID transition temperature is ksT} — 9.67fiaj z , corre- 
sponding to 0.269T C . 

To understand the significance of the 3D transition 
temperature T c we calculate the specific heat C, C = 
(dU/dT)^, where U denotes the internal energy [20j, for 
the ideal gas in the grandcanonical ensemble. The inset 
of Fig. n shows the specific heat C for N — 27 for three 
different aspect ratios, i.e., L — 1 (solid line), 10 (dot- 
ted line) and 100 (dashed line). Since the specific heat 



shows a peak, although broadened due to the finite size 
of the Bose gas, at T/T c w 1 for L = 1 and 10, and 
at T/T c ks 1.4 for L = 100, it is justified to speak of a 
3D transition temperature for Bose gases with as few as 
N = 21 atoms. In contrast, the transition to macroscopic 
occupation of the lowest energy state for the quasi-lD 
gas with L = 100, i.e., to "ID condensation", does not 
imprint a clear signature on the specific heat (see also 

Ref. mmi). 

Figure [5] shows the expectation value of \z\ in units of 



1.4 

1.2 

1.0 

0.8 

0.6 
3.0 
2.5 



JS 2.0 



N 

<a 

\ 

A 

N 

OB 

V 



1.5 
1.0 



A 

N 
CO 

XI 

<a 0.5 

v 

0.0 



6 
5 

4 

3 

2 

1 







N 

CO 

\ 

A 

N 

m 

ecl 

V 



(a) L=l * 
6° 

. o 



.0 



ooo 



5 < 



(b) L=10 a 



"5" 



A 

A 9 

a A 9*1.6 



>9o 



9 



o 1.4 

1.2 

1.0 
0.8 



6 



0.0 0.5 1.0 1.5 



(c) L=100 A 

A 

A g 



9 



9 



8 1.00 

n * 0.95 
o 

0.90 



Q.Q 0.5 l.Q 1.5 



0.5 1.0 1.5 
Temperature T/T c 



FIG. 2: PIMC expectation value of \z\ in units of a z as a 
function of T/T c for TV = 27 for (a) L = 1, (b) 10 and (c) 
100. Diamonds show the results for a/a z = 0, squares those 
for a/a z = 0.00433 and triangles those for a/a z = 0.0433. For 
L — 10 and L = 100, the insets show the expectation value of 
p in units of function of T/T c . Statistical uncertainties 

are smaller than the symbol size [2^1 . 

a z , calculated using the PIMC method, as a function of 
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the scaled temperature T/T c for (a) L = 1, (b) L = 10 
and (c) L = 100 for three scattering lengths; a — (dia- 
monds), a — 0.00433a z (squares) and a = 0.0433a z (tri- 
angles). At low T/T c , our expectation values of \z\ for 
a = (circles) approach the zero temperature value, i.e., 
(\z\) = 0.564a z . Since the energy of the transverse exci- 
tations increases with increasing L, the expectation value 
of \z\ for L = 100 approaches the zero-temperature value 
at a lower scaled temperature T/T c than that for L = 1. 
For repulsive interactions, i.e., a/a z — 0.00433 (squares) 
and a/a z = 0.0433 (triangles), the expectation value of 
\z\ increases compared to that of the non- interacting gas. 

The insets of Figs. EJb) and (c) show the expectation 
value of p in units of a p , where a p = (mw p ), for 
L = 10 and 100, respectively, as a function of T/T c (using 
the same symbols as in the main figure). At T = 0, the 
expectation value of p is 0.886a p for the non-interacting 
gas. Just as the expectation value of \z\, that of p de- 
pends strongly on the interaction strength a. The nearly 
constant expectation value of p for L = 100 (for a given 
value of a) at low T/T c indicates that the excitations in 
the transverse direction are frozen out for T < 0.4T C , and 
hence for T < T} D . 

We now turn to the calculation of the superfluid frac- 
tion (n s /n)n with respect to the axis n. For the non- 
interacting gas, the superfluid fraction with respect to, 
e.g., the z-axis, can be calculated from the thermal ex- 
pectation values of x 2 and y 2 0. To this aim, we con- 
sider a trapping geometry with oj y = lu x + Aw in the 
limit Alo -> (see Eq. (7) of Ref. 0). Dotted lines in 
Fig-IUshow the resulting superfluid fraction (n s /n)g, cal- 
culated in the grandcanonical ensemble, for (a) L = 1, 
(b) L = 10 and (c) L = 100 as a function of T/T c for 
N = 27 non- interacting bosons. 

Within the PIMC formulation, the superfluid fraction 
(n s /n)n can be calculated from the square of the pro- 
jected area A n . 1 Q| , where An = A-n and A denotes the 
area enclosed by the imaginary time paths @. Symbols 
in Fig. |2| show the superfluid fraction {n s /n)z calculated 
using the PIMC method for three different aspect ratios 
L] diamonds show our results for a = 0, squares those 
for a — 0.00433a z , and triangles those for a = 0.0433a z . 
For a = 0, the PIMC results for (n s /n)z (diamonds), cal- 
culated in the canonical ensemble, agree well with those 
calculated in the grandcanonical ensemble (dotted lines). 

The superfluid fraction (n s /n)i is essentially one at 
small scaled temperatures, and decreases gradually with 
increasing T/T c . For L = 1, {n s /n) z is about 0.05 for 
T/T c = 1. For the larger aspect ratios [see Figs. Efb) 
and (c)], in contrast, (n s /n)z is significantly larger at the 
transition temperature (about 0.2 for L = 10 and about 
0.65 for L = 100). When plotted, as done here, as a 
function olT/T c the superfluid fraction (n s /ri)g shows a 
very weak, if any, dependence on the interaction strength 
for all aspect ratios. 

The spherically symmetric system with L = 1 has no 
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FIG. 3: Superfluid fraction (n s /n)^ for N — 27 and (a) 
L = 1, (b) L — 10 and (c) L — 100: Diamonds show the 
PIMC results for a/a z = 0, squares those for a/a z = 0.00433 
and triangles those for a/a z = 0.0433 [23|: dotted lines show 
(n a /n)s for a = calculated in the grandcanonical ensem- 
ble. Solid lines show the superfluid fraction given by Eq. 
Dash-dotted lines show the condensate fraction No/N, and 
dot-dot-dot-dashed lines the fraction Nm/N. The inset of 
panel (a) shows {n s /n)± for N — 27 and T/T c — 0.2 for three 
different interaction strengths (using the same symbols as in 
the main figure) as a function of L on a logarithmic scale. 



preferred symmetry axis, implying (n s /n) z = (n s /n)x- 
For L > 1, however, the superfluid fraction (n s /n)$ is 
distinctly different from (n s /n) z . Geometric arguments 
imply that {n s /n)x approaches zero when the system 
reaches the quasi- ID regime. When exposed to a rota- 
tion about x, the atoms move with the external trap, thus 
implying 9" 9 = O^. The inset of Fig.^a) shows the su- 
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perfluid fraction (n s /n) & for TV = 27 and T/T c = 0.2 for 
three different scattering lengths (using the same sym- 
bols as in the main figure) as a function of the aspect 
ratio L. It is evident that (n s /n)x decreases rapidly with 
increasing L. 

To connect with earlier work, we consider an analytical 
expression for the superfluid fraction (n s /n) z , which has 
been derived using the semi-classical approximation for 
the non-interacting gas [?|]. For a trap geometry with 
ujy w tu x and us z = u> x /L, we generalize the treatment by 
Stringari to account for rotations about n = x. To 
additionally improve the accuracy for small N, we use 
the T c that accounts for finite-size effects (see above), 
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Here, A and B denote constants depending on the ge- 
ometry of the trap; A = B = 1 for n = z, and 
A = 1 - [(1 - L)/(l + L)] 2 and B = (L 2 + 1)/(L + 1) 
for n — x. Solid lines in Fig. [3] show the resulting ap- 
proximate superfluid fractions (n s /n)j (main figure) and 
{n s /n)x [inset of Fig. OJa)]. The agreement between 
Eq. (J2J and our numerical results for L = 1 is good. 
For L = 100, however, Eq. describes the superfluid 
fraction only qualitatively. In particular, Eq. J2J) clearly 
underestimates (n s /n) z for T/T c > 1. 

For comparison, dash-dotted lines in Fig. [21 show the 
condensate fraction Nq/N calculated in the grandcanon- 
ical ensemble for the non- interacting gas. For L = 1, the 
condensate fraction roughly agrees with the superfluid 
fraction. For L — 100, in contrast, Nq/N drops to zero at 
much lower temperatures than (n s /n) z . This shows that 
the condensate fraction and the superfluid fraction are 
distinctly different quantities for highly-elongated sys- 
tems. For comparison, dot-dot-dot-dashed lines show 
the fraction of atoms Nm/N in the lowest transverse 
mode, where Nm — X^fc-^oofc and N^k/N denotes the 
fraction of atoms in the state with i quanta in the x-, j 
quanta in the y- and k quanta in the z-direction. For the 
highly-elongated gas with L = 100, the fraction of atoms 
Nid/N is larger than the superfluid fraction (n s /n) z but 
shows a similar overall behavior. Calculations for the 
non-interacting gas for N = 1000 and L = 5000 (not 
shown) indicate similar behaviors. 

In sumary, this Letter describes microscopic calcula- 
tions for small Bose gases over a wide temperature range. 
We study the crossover from 3D to ID by changing the 
trapping frequency. Since all our calculations are per- 
formed in full 3iV-dimensional configuration space, the 
freezing of the radial motion at low temperatures emerges 
from our calculations; it is not an input. Specifically, 
we determine the temperature dependence of the quan- 
tum mechanical moment of inertia. This quantity has 
played a key role in the study of finite-size bosonic he- 



lium droplets over the past 10 years or so |2J| . Measure- 
ments of the quantum mechanical moment of inertia of 
an impurity embedded inside such a droplet have, e.g., 
shown unambigiously that bosonic helium clusters with 
as few as about 60 atoms are superfluid [25j . This paper 
shows that the effects of superfluidity are altered as the 
effective dimensionality of the trapped gas changes from 
3D to ID. The superfluid fraction (n s /n) z is enhanced as 
the dimensionality is reduced. In the quasi-lD regime, 
the superfluid response is distinctly different from the 
condensate fraction Nq/N and very roughly follows the 
fraction Nm/N of atoms in the lowest transverse mode. 
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